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Abstract 

After a brief review of various mappings of fermion pairs to bosons, we rig- 
orously derive a general approach. Following the methods of Marumori and 
Otsuka, Arima, and lachello, our approach begins with mapping states and 
constructs boson representations that preserve fermion matrix elements. In 
several cases these representations factor into finite, Hermitian boson images 
times a projection or norm operator that embodies the Pauli principle. We pay 
particular attention to truncated boson spaces, and describe general methods 
for constructing Hermitian and approximately finite boson image Hamiltoni- 
ans, including effective operator theory to account for excluded states. This 
method is akin to that of Otsuka, Arima, and lachello introduced in connec- 
tion with the Interacting Boson Model, but is more rigorous, general, and 
systematic. 
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I. INTRODUCTION 



The original attempt at bosonization of the nuclear many-fermion system was motivated 
by collective particle-hole modes in nuclei Since that time the interacting boson model 
(IBM) has been phenomenologically very successful in explaining low energy nuclear 
spectroscopy for heavy nuclei. The bosons in this model are thought to represent monopole 
(J=0), quadrupole (J=2), and sometimes hexadecapole (J=4) correlated pairs of valence 
nucleons in the shell model. The IBM Hamiltonian is Hermitian, usually has at most two- 
boson interactions, and conserves boson number, reflecting the particle-particle, rather than 
particle-hole, nature of the underlying fermion pairs. While one can numerically diagonalize 
the general IBM Hamiltonian, one of the strengths of IBM is the existence of algebraic 
limits corresponding to the subgroups SU(3), U(5), or 0(6), with analytic expressions for 
excitation bands and transition strengths, which encompass an enormous amount of nuclear 
data. 

The microscopic reasons for the success of such a simple model are elusive. Otsuka, 
Arima, and lachello, along with Talmi, have used a mapping of the shell model Hamiltonian 
to the IBM Hamiltonian based on the seniority model 0, but these attempts have not 
done well for well-deformed nuclei 0. For this reason we have revisited boson mappings to 
see if we can understand the success of the IBM starting from the shell model. 

We will begin by sketching out various historic approaches to boson mappings 
We then specificially follow Marumori and Otsuka et al. (OAI) in our mapping 
procedure which maps fermion states into boson states and construct boson operators that 
reproduce fermion matrix elements. We give the boson representation of the Hamiltonian 
and review the result that, in the full Boson Fock space, it factorizes into a boson image, 
which is the same as the Belyaev-Zelevinskii Hamiltonian [Q-in fact in the full space all 
mappings yield the same results- times a normalization operator which projects out the 
spurious states. However, since our goal is to understand the IBM, which only deals with a 
few of the enormous degrees of freedom of the shell model, we go on to discuss boson images 
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in truncated spaces. This, we shall see, gives rigorous insight into the OAI mapping and 
shows how to systematically extend it. 



II. A BRIEF HISTORY OF BOSON MAPPINGS 

The fundamental goal is to solve the many-fermion Schrodinger equation 

i^l^A) =^a|^a) (1) 

and find transition matrix elements between eigenstates, t\x' = (^^\ 
many-body (shell-model) basis states one often uses Slater determinants, antisymmetrized 
products of single-fermion wavefunctions which we can write using Fock creation operators: 
■ ■ ■ |0) for n fermions. For an even number of fermions one can instead construct states 
from products of fermion pairs, 

i*^)=n4jo); (2) 

m=l 

if the number of fermion is fixed at n then m runs from 1 to = n/2. As the fermion Fock 
space may be so large as to make direct solution intractable, the idea of a boson mapping 
is to replace the fermion operators with boson operators, using only a minimal number of 
boson degrees of freedom, that approximate the spectrum and transition matrix elements of 
the original fermion problem. There are two approaches to boson mappings which we now 
review. 

The first approach, epitomized in nuclear physics by Belyaev and Zelevinskii (BZ) is 
to map fermion operators to boson operators so as to preserve the original algebra. Specif- 
ically, consider a space with 2f2 single-fermion states; gnify fermion creation and 
annihilation operators. The set of all bilinear fermion operators, form the 
Lie algebra of SO(4f2), as embodied by the commutation relations 

[aittj^akai] = (3) 
= SiiSjk + Sikalaj + 6jialai - {i ^ j) (4) 
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T ). For fermion 



= Sjkttiai - (i ^ j) 
= Sjkajai — 6iial.aj 



(5) 
(6) 



At this point it is convenient to introduce collective fermion pair operators 



We always choose the Q(2Q — 1) matrices Ap to be antisymmetric so as to preserve the 
underlying fermion statistics, thus eliminating the need later on to distinguish between 
'ideal' and 'physical' bosons. We also assume the following normalization and completeness 
relations for the matrices: 



tr AqA^ — ^af3] 



(8) 
(9) 



■tj ~ -J - 2 

Generic one- and two-body fermion operators we represent by T = J^ijTijalaj, V = 
S;ui/ (a*I ^ 1^) ^ji^i'j where Tij = {i\ T from such operators one can construct a fermion 
Hamiltonian H. Now one has the following commutation relations: 



= [iLiy =0; (10) 
A^, Al] = - 2 ^ (aJjA^) . . ajaj; (11) 

'i^, r] = 2 ^tr (A^TA^ Ap (12) 

% f^l = E (T1T2 - T2T1),, ajaj (13) 

ij 

The method of Belyaev and Zelevinskii is to find boson images of the bifermion operators, 



i^i)B = &i + E ^f%bib, + E xf^Hibibib.b, +... 

a/37 ap'ySe 

{f)B = Ey''%bp + E y^'^%b\b0bs + ... 



(14) 

(15) 
(16) 



where 6^, b^p are boson creation and annihilation operators, [ba, = ^a/j, with the coefficients 
x,y chosen so that the images {A^^)b, (^i/)b, {T)b have the same commutation relations as 

4 



in ([T0|)-([r3D. Because the algebra is exactly matched, if one builds boson states in exact 
analogy to the fermion states then the full boson Fock space is not spanned and one does 
not have nonphysical or spurious states. 

In the full boson Fock space, that is, no truncation of the boson degrees of free- 
dom, the image of one body operators is finite and given quite simply by (T) b = 
2J2^ptT{AaTAl)blbp. Since in the full space any fermion Hamiltonian can be written 
in terms of one-body operators, the boson image of a finite fermion Hamiltonian will be 
finite in the full space. The states that one must use then are built from the boson repre- 
sentations of the fermion pairs given in (|l^) which will not just be products of bosons but 
will include exchange terms. For example, for two bosons and using (p!4D, 

AiAUo) (bibu xi-^ibi) \o) . (17) 



These exchange terms are due to the antisymmetry. We shall take care of such exchange 
effects by introducing a norm operator in the boson space. 

For truncated spaces, however, blind application of the boson representations of the pair 
operators (|l^ will produce states outside the truncated space. Therefore, if these states 
are omitted, then the Hamiltonian needs to be renormalized to account for these omissions, 
which in general leads to an infinite BZ Hamiltonian. Marshalek points out that there 



exist mappings that are both finite and Hermitian, but these in general require projection 
operators to eliminate spurious states. We will regain this result later on in this paper. 
A variant of Belyaev-Zelevinskii that also preserves commutation relations is the Dyson 



mapping |1 1|| : 



i„ ^ (18) 
4 _ 6^ _ 2 5: tr (A,At A,At )6l6t 6. (19) 

f ^ 2tr Y.{A^TAl)b%. (20) 

The operators are then clearly finite; on the other hand they are just as clearly non- 
Hermitian. From a computational viewpoint non- Hermit icity is only a minor barrier, but 



it is an obstacle to an understanding of the microscopic origin of Hermitian IBM Hamilto- 
nians. Furthermore, unless the truncated pairs constitute a subalgebra, under truncation, 
renormalization will produce an infinite expansion, although it may be possible that this 
expansion may be more convergent than the BZ expansion. 

The second major approach, pioneered by Marumori 0, is to map fermion states and 
construct boson representation operators that preserve matrix elements. The original work 
of Marumori, however, focused on particle-hole excitations and so the number of pairs and 
consequently bosons were not fixed. However, this method can be applied to particle-particle 
pairs as well. 

Marumori constructs the norm matrix 

Ar„^ = (^„|^^) (21) 

and then the Usui operator 

f/= E \'^p)i^r/f{^o.\ (22) 

where bosons states are constructed in strict analogy to the fermion states, 

N 

l*)= n^LlO)- (23) 

m=l 

Then the Marumori expansion of any fermion operator is 

Ob = UOpUK (24) 

We shall show that, in the full space, these boson representation operators factor into a finite 
boson image times a norm operator. Furthermore, the boson image of the Hamiltonian is 
the same as the BZ image, and hence, the two methods are equivalent in the full space. 

Otsuka, Arima, and lachello (OAl), along with Talmi investigated the microscopic 
origins of the Interacting Boson Model through boson mappings. Although they also mapped 
states, they differed from Marumori in some key details. First of all, they built states built 
on a fixed number of particle-particle, not particle-hole, pairs. In addition, the space was 
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truncated to include only one monopole (J'^ = 0"*") and quadrupole {J'^ = 2"*") pair. These 
states were orthogonalized based on seniority. That is, they constructed, for 2N fermions, 
low-seniority basis states of S and D fermion pairs, S^^'^'^D"-'^'^, and then orthonormalized 
the states such that the zero- seniority state is mapped to itself, and states of higher seniority 
V were orthogonalized against states of lower seniority, 

\v) I "t;") = \v) + \v-2) + \v-A) + ... (25) 



Then OAI calculated the matrix elements ( "^^""dD^'d" \Hf\ for nd^n'^ = 

0,1,2 and obtained the coefficients for their one plus two-boson Hamiltonian. These coeffi- 
cients have an implicit A^-dependence (and for large and arbitrary systems such matrix 
elements are not trivial to calculate, especially in analytic form!) and thus a many-body 
dependence. At first sight this is not entirely unreasonable as it is well known the IBM 
parameters change substantially as a function of the number of bosons, even within a major 
'shell'. Nonetheless the OAI mapping has three drawbacks. The first is that it is not clear 
how to systematically calculate many-body contributions beyond that contained in the OAI 
prescription, whereas the method we shall describe is fully and rigorously systematic. The 
second is that the OAI prescription can induce many-body effects where none are needed. 
This point will be illustrated in section [V Q . Thirdly, only the = 0, 1, 2 space is exactly 



mapped, but very deformed systems will involve large n^. In fact, for an axial rotor limit, 
the average number of d-bosons in the ground state band is 2/3 the total number of bosons. 
Correcting this by also mapping matrix elements with > 2 will involve many-body terms. 
As an alternative to OAI, Skouras, van Isacker, and Nagarajan |]12[ proposed a "demo- 



cratic" mapping where the orthogonalization is based on eigenvectors of the norm matrix 
rather than seniority. 

In what follows we present a unified state-mapping method and obtain four strong results. 
First we derive matrix elements of the fermion operators in the pair basis (2). Second, we 
give general expressions for fermion matrix elements via boson representations. Third, we 
show how in several cases one can have exact, finite, and Hermitian boson images of fermion 



operators. Finally, we show how to extend the OAI and democratic mappings in a systematic 
and rigorous fashion, and illustrate how the choice of orthogonalization can affect the many- 
body dependence of the boson images. 



III. MATRIX ELEMENTS OF FERMION-PAIR STATES 



The starting point of any state-mapping method is the calculation of matrix elements of 
fermion operators between states constructed from fermion pairs of the form (^, including 



the overlap: (^al^/j), 



H 



T 



and so on. These matrix elements are 



much more difficult to compute than the corresponding matrix elements between Slater 
determinants. As we shall show, however, full and careful attention paid to the problem of 
calculation matrix elements can yield powerful results. Silvestre-Brac and Piepenbring [jl3 



laboriously using commutation relations, derived a Wick theorem for fermion pairs. Rowe, 



Song and Chen WM using 'vector coherent states' (we would say fermion-pair coherent states) 



found matrix elements between pair-condensate wavef unctions, states of the form [A^j |0) . 



Using a theorem by Lang et al. |]T5[, we have generalized the method of Rowe, Song and 
Chen and recovered (actually discovered independently) the expressions of Silvestre-Brac 
and Piepenbring. 

To derive the generalized Wick's theorem, we use the following theorem (its proof, which 
requires fermion coherent states and integration over Grassmann variables, is found in Ap- 



pendix A of [|15|] ): 

Let U be an operator of the form 



U = exp (h{n)) ■ ■ - exp {h{2)) exp (h{l)) 



(26) 



where each h{t) is of the form 



h{t) = \^T{t)ijaldj + B{t)ijaiaj + A*{t)jialal 



(27) 



Again, the matrices A, B are antisymmetric (T is not, in general), and are of dimension 
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2Q X 2Q. We introduce the 4Q x 4Q matrix representation of U in the basis of the single- 
particle Fock operators aj,aj. Then: 

Un Ui2\ n / T{t) 2At(t)\ 

= nexp ^ . (28) 

U21 V22J t=i V2B(t) -T^(t)/ 

The conclusion of the theorem is the vacuum expectation value of If : 

(0|f/|0) = v/det (U22) exp (^EtrT(t)) . (29) 



Appropriate derivatives of (p9D bring down pair creation and annihilation operators, 
leading to the desired matrix elements. We do this in some detail for the overlap matrix 
elements {9a\ ^/s) for the states given in (2); then we simply give the results for matrix 
elements of one- and two-body operators which are found in the same way. If one begins 
with 

N N 

[/ = n exp {e^X.) n exp (e/3,4j ' (30) 

i=l i=j 

where Aa, At are defined in (7), then the overlap is 



(31) 

alle's=0 



In applying the theorem the matrix algebra is straightforward; keeping only terms linear 
in ea- etc. one finds that for the overlap T(t) = 0, and 

U22 = l + 25:e„,^e^^A,,Aj^.. (32) 

ij 

Next, insert U22 into (pQ]), use det U22 = exptrlnU22 and expand the logarithm to arrive 
at the generating function 

(0| U |0) = exp (f: ^^^^tr e„.e^ A« A^ . ) ] . (33) 
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Finally, one applies (|3lD. 

The result is constructed from objects of the form tr ( A^^ Ajj^ Aq,2 Ajj^ . . . Aa^A^^], which 



we term a k-contraction (related to the contractions of k phonons in [|l^]) and represent with 



the notation Ck[(yi, Pi, a2, P2', ■ ■ ■ c^k, Pk] (from the order it is clear that the a's represent 
final states and the /3's initial states). Then \E'^) is a finite sum of products of k- 
contractions; each term in the sum is a product of mi 1-contractions, 1712 2-contractions, 
. . . , and rriN A^-contractions, with the {rrii} taking all possible values under the constraint 
that J2k kmk = N. For example, there will be terms consisting of 1-contractions, terms 
consisting of — 2 1-contractions and one 2-contraction, and so on, up to terms consisting 
of 1 A^-contraction. The full expression is 

N N 
(^,|^;3)=EE^°(«l'«2,...,«£;A,...,/3£)((Ar-^)!)-l n ClK,/?fc] (34) 
P £=0 k=e+i 

where V means the sum is over all permutations of {aj} and {Pj}- The coefficients are 

given by 

^'i = J2 Zi{m2,...,Tni) 

m2,...,Tni 

X C2[ai,Pi; a2,P2] C2[a3, P^; a^, P^] . . Pm,] (35) 

V ' 

1712 

X C3[a[, P[; a'2, P2; a'^, P'sl ■ ■ 

^ v ' 

The organization of (|34D is such as to ease the interpretation in bosons in the next section. 
The coefficients Zi{m2, ■ ■ ■ ,me) are found through the expansion of (p^): 

z,{m2,...,m,) = i-2Y'^^l[^^, (36) 

where i = J2k>i ^"^fc and M = J2k>i^k, and Wq = Zq = 1 (hence the normalization chosen 
above). Note that we have set mi = A^ — £ and that the w'^ have no explicit dependent on 
A^. 

To illustrate, consider a single j-shell fermion space, with fl = j ' + 1, and further consider 
only the fermion pair with total J = 0: Aq = (2)^^/^ 
agonal, (aJ) ^ = (2fi)"^/^(-l)-''~™'5m Then the ^-contraction Ck = tr (AqAJ,)" = 
{2fl)~''~^^. Following the above definitions, 

^2(1) = -1, ^3(0, 1) = 4/3, z,{2, 0, 0) = 1/2, 
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1] a] 



the matrix An is antidi 



Z4(0, 0, 1) = -2, Z5(l, 1, 0, 0) = -4/3, 25(0, 0, 0, 1) = 16/5, etc. 



and 



1 1 1 1 







1 



1 



(37) 



Then from equation (p4D, including the sum over all permuations which gives a factor A^!^, 
one obtains for example the norm of the state with 5 J = pairs: 







n\ r, /1 10 35 50 24- 



1-N/n {N 



which agrees with the general result from the commutation relation 
is the fermion number operator), that is, the norm is N\Q\/{Q — N)\Q^^^. 

In the same way as for the overlap one can derive the matrix elements for one- and 
two-body operators O^'^ (and for general n-body operator, if so desired): 



N 



V i=l 



e N 
-^E*fc'[^'''i«i'---'/^fcK-fc n (38) 

fc=l m=e+i 



For a one-body operator T we define Aq,T = A^m] then 



wl{T;ai,a2,...,Pk) = 2(-2)^ ^Ck[ai{T), Pi, a2, P2; ■ ■ ■ ; Oik, Pk\- 



(39) 



For a two-body operator V 



Wfc(F;ai,a2,...,/3fc) = 



fii/ 



1=1 



-2((5;,i - l)cfc+i[/i,/3i;a2,/32; • • • ; a^. A; "i, a^+i, A+i; • • • ;afc,/3fc]) . (40) 



Returning to our illustration, first consider the the number operator A^, which is represented 
by just the unit matrix; in this case, again restricting ourselves to J = pairs, wl = 
2/(— r2)^^^. Similarly for the pairing interaction V = AqAq, 



r,2 



W 



-n 



k-l 



k 1- 



The reader is invited to check that these coefficients reproduce the correct matrix elements. 
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Therefore, given any two states constructed from fermion pairs and the matrices repre- 
senting those constituent pairs, the above formulas give exactly the overlap and the matrix 
element for one- and two-body operators. For a pair-condensate wavefunction, the matrix 



elements can be found quickly through recursion [p!3| , p |. 

Throughout this paper we will discuss the effect of truncating the boson Fock space, that 
is taking a restricted number of boson species from which to construct states and operators, 
on boson mappings. But which boson species should we keep? In regards to this question we 
merely wish to comment that Rowe, Song and Chen give a variational principle which 



seems useful in this regard, and is probably related in some approximation to the Hartree- 
Fock-Bogoliubov states that Otsuka and Yoshinaga use in their mapping of deformed 
nuclei. 

IV. BOSON REPRESENTATIONS OF FERMION MATRIX ELEMENTS 

We now want to translate the fermion matrix elements into boson space. We take the 
simple mapping of fermion states into boson states 

N 

\^p)^\^p) = n^Lio)' (41) 

m=l 

where the 6^ are boson creation operators. We construct boson operators that preserve 
matrix elements, introducing boson operators Tb, Vb, and most importantly the norm 
operator jifs such that ($q|7b|$/3) = (^a|T|^/3), = (^olV^I*/?)- and 

(<I'q,|A/b|$/3) = 1^/3)- We term 7b,Vb the boson representations of the fermion op- 
erators T, V. The boson basis 1$^) is an orthogonal basis. The fermion norm operator in 
the boson space will be given by 

^B = T.^-^0^bl..bij,,...b,, (42) 

where the coefficients are given by (P^. Because the matrices Aq,, are orthogonal (see 
(8)), the one-contraction is simply Ci [ak,Pk] = ^a^fik- Using the fact that b'^Jia = the 
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number operator, we find the 'linked-cluster' (a la Kishimoto and Tamura p!7|JT8[] although 
with differences) expansion of the representations to be of the form 

oo It 
1=2 {a,T} i=\ 3 = 1 

and similarly for Vb,7b- In the norm operator the £-body terms express the fact that 
the fermion-pair operators do not have exactly bosonic commutation relations, and act to 
enforce the Pauli principle. 



In the example of a single j-shell given in the previous section, using the coefficients (pTf) 
and with the mapping ^ s^, the purely s-boson part of the norm is 

A^^ = 1 - ^shhs + ^.sUhUss + - ^) sUhhUsss + ... (44) 

which again yields the correct matrix elements. 

The norm operator can be conveniently and compactly expressed in terms of bosons by 
using the fermion norm matrix (^) in terms of the fermion generating function (^). Taking 
derivatives with respect to the e^j- in ([31| ) is like contracting 6^. with a b^-, and with respect 
to e^. is like contracting a ba^ with a fej^^. Hence the norm matrix is just the generating 
function ( p3D with — > and — * bjs.. However, since these bosons do not commute 
with one another, we must take the normal order: 

J^B=:exp(^-\p^t^C,y. (45) 

where the colons ':' refer to normal-ordering of the boson operators, and Ck = 2:trP'^: is 
the /cth-order Casimir of SU(2fi), with P = J^b^brA^Al (the trace is over the matrices 
and not the boson Fock space). This norm operator, which takes into account the exchange 
terms in the BZ expansion of a fermion pair given in (0), is found in Ref. [jl9 . 



Similarly — and this is a new result we have not seen elsewhere in the literature — the 
representations Tb, Vb can also be written in compact form: 

fB = 2j2: tr [A^TAt gI blbrJ^B: (46) 
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= E 1^1 ^) E- [A.At g] tr [A.At G 

+ 4tr [A^AtpGA^AtG] } hlKMB: , (47) 

where G = (1 + 2P)^^. These compact forms are useful for formal manipulation. Further- 
more they have the powerful property of exactly expressing the fermion matrix elements 
under any truncation, a fact not previously appreciated in the literature even for the norm 
operator [0. By this we mean the following: suppose we truncate our fermion Fock space 
to states constructed from a restricted set of pairs {a}. Such a truncation need not corre- 
spond to any subalgebra. Then the representations in the corresponding truncated boson 
space, which still exactly reproduce the fermion matrix elements and which we denote by 
[Mb\t etc., are the same as those given above, retaining only the 'allowed' bosons with 
unrenormalized coefficients. For example 

WB]r =■■ exp ( -- E ^ \Ck\ ^ ) : (48) 



where 

[C,]t = 2:tr ([P]^)'=: , [P]^ = E4&.A^At. (49) 

CTT 

This invariance of the coefficients under truncation will not hold true for the boson images 
introduced below. 

With the boson representations of fermion operators in hand, one can express the fermion 
Schrodinger equation ([^) with = T + as a generalized boson eigenvalue equation, 

HB\<^x)=ExAfB\^x). (50) 

Here Ti.B is the boson representation of the fermion Hamiltonian. Every physical fermion 
eigenstate in (|l|) has a corresponding eigenstate, with the same eigenvalue, in (^0]). Because 
the space of states constructed from pairs of fermions is overcomplete, there also exist 
spurious boson states that do not correspond to unique physical fermion states. These 
spurious states will have zero eigenvalues and so can be identified. 
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V. BOSON IMAGES 



In general the boson representations given in (^5]), (^61) and (^Tj) do not have good 
convergence properties, so that simple termination of the series such as (^31) in £-body terms 
is impossible and use of the generalized eigenvalue equation (|50|), as written, is problematic. 
Instead we "divide out" the norm operator to obtain the boson image, i.e. schematically, 

h ~ "Tiij/Mj." (51) 

That this is reasonable is suggested by the explicit forms of ( ^Bj ) and (^). The hope of 
course is that h is finite or nearly so, so that a 1+2-body fermion Hamiltonian is mapped 
to an image 

h ~ 0ib^b + 02b^b'^bb + O^b^b'^b'^bbb + O^b'^b^b^^bbbb + ... (52) 

with the £-body terms, £ > 2, zero or greatly suppressed. We now discuss how to "divide 
out" the norm. 



A. Exact results: Full Space 

It turns out that for a number of cases the image of the Hamiltonian is exactly finite. 
In particular, for the full boson Fock space the representations factor in a simple way: 
Tb = N'bTb = TbN'b and Vb = J^syB = Vb-^b, where the factored operators Te, Vb, 
which we term the boson images of T, V, commute with the norm operator and have simple 
forms: 

fB = 2^tr [A,TAl)bibr, (53) 

VB = Y.if^\yw) 

The proof of the factorization and commutation requires use of the identities 



bib. + 2EEtr (A^At A.,At A,A|,) 6^6^6.6.. 



(54) 
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2 tr(QA|^)tr(A„R) = tr(QR) - tr(Q^R), (55) 

a 

2 tr(QAj,RA,) = tr(Q)tr(R) - tr(Q^R) (56) 

a 

(Q^ is the transpose of Q) which in turn are proved using the completeness relation (^. 
The image Hamiltonian Hb = Tb + Vb is the one determined by BZ if one decomposes the 
Hamiltonian into multipole-multipole form and then maps these multipole operators. As 
discussed earlier, these BZ multipole operators are finite in the full space. This result, and 



its relation to other mappings, was noted by Marshalek @,10 

Thus any boson representation of a Hamiltonian factorizes: Hb = MbHb in the full 
space. Since the norm operator is a function of the SU(2n) Casimir operators it commutes 
with the boson images of fermion operators, and one can simultaneously diagonalize both 
Hb and A/^. Then Eqn. (pDj ) becomes 

Hb\^x) = E'^\^x). (57) 

where E"^ = E\ for the physical states, but E'^^ for the spurious states is no longer necessarily 
zero. The boson Hamiltonian Hb is by construction Hermitian and, if one starts with at 
most only two-body interactions between fermions, has at most two-body boson interactions. 
All physical eigenstates of the original fermion Hamiltonian will have counterparts in (pT]). 
It should be clear that transition amplitudes between physical eigenstates will be preserved. 
Spurious states will also exist but, since the norm operator Mb commutes with the boson 
image Hamiltonian Hb, the physical eigenstates and the spurious states will not admix. 
Also the spurious states can be identified because, while they will no longer have zero 
energy eigenvalues, they have eigenvalue zero with respect to the norm operator. 



B. Exact Results: Truncated space 

The boson Schrodinger equation (|57D, though finite, is not of much use as the boson 
Fock space is much larger than the original fermion Fock space, and we still must truncate 
the boson Fock space. Although the representations remain exact under truncation, the 

16 



T 



Me 



Again 



factorization into the image does not persist in general: Tii 
for example consider the pairing interaction, even in multiple j-shells, and a truncation 
to just s-bosons; then [Hb]t = s^s + ^s'^s'^ss whereas the appropriate image is actually 



h 



-s^s^ss. This was recognized by Marshalek [|T^. (An alternate formulation 



does not require the complete Fock space, but mixes physical and spurious states and so 
always requires a projection operator.) 

If the truncated set {a} represents a closed subalgebra, that is, if the truncated set of 
fermion pairs are closed under double commutations: 



~ ^a,(3 ~ Tap, 



(58) 
(59) 



then a factorization [|T^ 



(60) 



does exist, with /i^ at most two-body, but not necessarily Hermitian: 



err 

(4Ap)^ = 6t [ftp - E r?,.,4M.') 

\ arr' / 



(61) 
(62) 



We term this a generalized Dyson image PJ9|JTT[]. Although (physical 



h 



D 



spurious 



0, if hi) is non-Hermitian, that is A/s 



h 



D 



7^ 0, then unfortunately 



spurious 



D 



physical^ 7^ 0. In the full space, of course, all definitions of boson images 
coincide and yield the same result. 

We have found conditions under which ho is additionally Hermitian and commutes with 
the truncated norm operator A/^ • Consider a partition of the single fermion states labeled 



T 



by i 



where the dimension of each subspace is 2Qa, 2fic so that Q = 2Qa^c- We 



denote the amplitudes for the truncated space as and assume they can be factored, 



fA, 



fKt), 



{Al)i^j^, with KtK = KK"^ = ^ and = {-IfK, where p = 
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(symmetric) or j9 = 1 (antisymmetric). Furthermore we assume the completeness relation 
(y), which was crucial for proving that Hb = J^bHb 0, is valid for the truncated space; 
i.e., 

The norm operator in the truncated space then becomes 

N.--pg(^)""i*^(p')^. («) 

where P = Yl,af ^s-^f A^-A^ so that [P]j, = [i^f^ P. In this case the boson image of a 
one-body operator is the truncation of the boson image in the full space, 



B 



Ml 



B 



T, 



B 



(65) 



'-B 



2^tr (A^TAt) 46,. 



(66) 



The representation of a two-body interaction can be factored into a boson image times the 
truncated norm. 



B 



however, vb, while finite (1+2-body), is not simply related to 
one-body operators. If one writes 

VD = Y.{^\v\f)hlh,+ (aa>|ff')44,M.', 

fjT aa'ff' 

then matrix elements of the two-boson interaction are 



B 



(67) 



as is the case for 



(68) 



(crcr'|f |rr') = 
-tr,{tre(A^A^A,A^,)tr,(A^,At; 



5^,(2^],-(-l)p)(fi, + (-l)f)' 

+2fiJtr,(A^A;^A^,At A,A|,) - tr,(A^AiA^,At ApKA^O] 

-^]„(2^], + (-l)'')tr,(A,KtAt A^,A|,)5^,m}, 



(69) 



where tia, tic are traces only in the a- and c- spaces, respectively. 



Upon inspection one sees the image ( p9| ) is not constrained to be Hermitian. Consider 
the additional condition between the matrix elements of the interaction: 



i;(f|Vk>E(A„k.,,.,,(At 

lJ..y la,3a '--l' 



(70) 



?/ 7' 



where the factor A^^^ = VLa{'^^a + (—1)^) is the number of pairs in the excluded subspace. 



While condition ([70| ) looks complicated there are interactions that satisfy it; for example, 
two-body interactions constructed from one-body operators V = T^pT^'p' where T^p = 



Al„Af^ . When (|70D is satisfied then vd is Hermitian and although vd ^ 
simply related: 

(T,r 

+2^^5:(/.|F|z.) Y: tr (A,AiA,,At A.At,) 44,M.' 



V, 



B 



they are 



V, 



B 



(71) 

by a factor which ranges 



with J'q^ = AQ^^/Na renormalizing the two-boson part of 
from unity (full space) to 2 for a very small subspace. Not all interactions satisfy ([TOD; 
for example, the pairing interaction never does except in the full space. For the pairing 



mteraction (/i 



ypairing 



S^l,oSu,oG, and AqAJ = and the image ( |67D ■yPj^"'™^ i^ecomes 



(remembering Q = 2Qa^ 

gUoIi- ^ 



N + ^ + ^]+ E tr (A,At AoAt,)46SM.' 



(72) 



where N is the total number of bosons, N = J2f blb^, and Nq = ^o^o- The second term in ( |72D 
is not Hermitian but can be transformed away by a similarity transformation leaving 
the first term as a finite Hermitian image which gives the correct eigenvalues for all N. 

The SO (8) and Sp(6) models |^ belong to a class of models which have a subspace 
for which (|63|) is valid and interactions which satisfy (|70|). In these models the shell model 
orbitals have a definite angular momentum j and are partitioned into a pseudo orbital angular 
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momentum k and pseudospin i, j = k + i. The amplitudes are then given as products of 
Clebsch- Gordon coefficients, 

(aI),, = {kmi,kmj\KaMa) (i/ij, /Xq), (73) 

where K and I are the total pseudo orbital angular momentum and pseudospin respectively 
of the pair of nucleons. For the SO (8) model i = | and one considers the subspace of 
pairs with K = {p = 0), {Al,)ij = ^'^^^^^^ \ i fii,i fij\Ia fia)i in the Sp(6) model k = 1 
and one considers the subspace with / = (p = 1), = ^^^^^^^ \ kmi, kmj\KaMa). 

The complicated conditions ( [70D hold true for important cases, such as the quadrupole- 
quadrupole and other multipole-multipole interactions in the S0(8) and Sp(6) models (that 
is, interactions of the generic form ■ in the notation of p2[) where therefore have 
Hermitian Dyson images. Not all interactions in these models have Hermitian Dyson images. 
For example, pairing in any model (see ([72|)) and, in the SO (8) model, the combination 
^7 = ^t^ + lp2 . p2^ ^j^g^g g ^ v^Aj=o, Vl = Ak + 2, which is the S0(7) hmit. It so 
happens that these particular cases nonetheless can be brought into finite, Hermitian form 
as discussed in the next section. 



C. Approximate or numerical images 



The most general image Hamiltonian one can define is 



h = U 



B 



-1/2 



T 



n 



B 



-1/2 



(74) 



which is manifestly Hermitian for any truncation scheme and any interaction, with U a 
unitary operator. (Because the norm is a singular operator it cannot be inverted. Instead 



B 



-1/2 . 



is calculated from the norm only in the physical subspace, with the zero eigenvalues 



which annihilate the spurious states retained. Then h does not mix physical and spurious 
states.) If W = 1, this is the democratic mapping [p!2|. Again, for the full space N'b,Hb =0 
and hence h = Jid = Hb- 
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This prescription is, we argue, useful for a practical derivation of boson image Hamil- 

tonians. Ignoring for the moment the unitary transformation W, consider the expansion 

r " 1 r ~ 1 ^1/2 

(^) of h. The operators 'Hb ^ and Mb ^ have similar expansions, and by multiplying 
out ( |7^ one sees immediately that the coefficient 6£ depends only on up to £-body terms 

r ~ 1-1/2 

in [HB]rp and A/^ ^ , derived from 2£-fermion matrix elements which are tractable for i 
small. Ideally h would have at most two-body terms, and our success in finding finite images 
in the previous section gives us hope that the high-order many-body terms may be small; at 
any rate the convergence can be calculated and checked term-by-term. Specifically, consider 
the convergence of the series (^) as a function of i. A rough estimate is that, for an A^- 
boson Fock space, one can truncate to the £-body terms if for i' > i, Opj is sufficiently small 
compared to 6^ x (A^ — £')!/( A^ — £)!; the strictest condition is to require O^i <^ Qijijl' — t)\. 

Although we have given in section IV analytic expressions for the coefficients of the 
boson representations, in practice one only needs the fermion matrix elements, given in 
section III, for the norm and the Hamiltonian or other operators. The coefficients of the 
image h are then found by numerical induction. Suppose one has the coefficients of h up 

-1/2 



AT; 



-note that in calculating 



to ^-boson terms. One takes the matrix elements of 

- ~ 1 -1/2 

Mb ^ , one first truncates and then calculates the inverse-square- root; the two operations 
do not commmute! — and of the Hamiltonian or transition operator in the 2{i + 1) fermion 
space, and multiply out those matrices as in (|7^ , yielding the matrix elements of h in the 
i + 1-boson space. The coefficients 9i+i of h in (|5^) are then uniquely determined, up to the 
freedom embodied in the unitary transform U. 



The Hermitian image h, defined in ([7^), is related to the Dyson image ho, defined in 
(|60|) , by a similarity transformation S =U Mb ^ , 

h = ShDS-\ (75) 

The similarity transformation S orthogonalizes the fermion states \^a) inasmuch 
{S-^^MbS-'^ = 1 in the physical space (and = in the spurious space). This is akin 
to Gram-Schmidt orthogonalization and the freedom to choose W, and iS, corresponds to 
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the freedom one has in ordering the vectors in the Gram-Schmidt procedure. The OAI and 
democratic mappings are just two particular choices out of many; the former orders the 
states by seniority whereas the latter takes U = 1. 

We can use this freedom in the choice of U to our advantage, which we illustrate in the 
SO (8) model [^], truncating to the space oi K = pairs (see eqn. ([75D), for which one 
has only an / = J = pair (mapped to an s boson) and a quintuplet of / = J = 2 pairs 
(mapped to dm bosons, m = —2 ... 2). To second order, the norm is 



B 



2Q 



1 

2fi 



L=0,2,4, 



d(S)d 



V5 r ^ 



d®d 



+ 



S ® d^ 



ss \ — 



ridUs 



(76) 



We pay particular attention to three interactions which correspond to algebraic limits: the 
pure pairing interaction Yv^^™& = S'^S = QAqAq, the quadrupole-quadrupole interaction 
yQQ ^— p2 . p2 -j-^ notation of |^), which can be written in terms of SO (6) Casimir 
operators, and the linear combination of pairing and quadrupole V'^ = \/p^'"'^s + i which 
can be written in terms of SO (7) Casimirs. As discussed in the previous section, the Dyson 
image of \/*3<9 jg Hermitian and finite, and hence hi) = h with U = 1: 



d®d 



+ 



d'^ ® d'^ ^ ssj + Sndfig 



(77) 



The Dyson images of the pairing and SO (7) interactions are finite but non-Hermitian: 

j^^pairing^^ = Qs^ S - S^^SS - V^S^^ [^^ ® ^] q " ^^^fls, 



(y'^^^ = {n + 5)sh + d^ ■ d-shUs + (g) S 



SS. 



(78) 
(79) 



We arrived at these images by computing the fermion matrix elements using the methods of 



2^1 and writing the norm and the representations of the interactions as matrices, which in 



second order are at most 2x2, and then directly multiplied the matrices JV ^Ti. We went to 
third order to confirm the Dyson images are finite. Using these same matrices we could also 
calculate M-^I'^UM-^I'^ , taking W = 1, which yields again (^ for V^^\ for 1/pai"°g and 
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the images are then Hermitian but with nonzero third-order, and presumably higher-order, 
terms. 

We then found Ws 7^ 1 for both the pairing and SO (7) cases (but not the same U) such 
that their respective Hermitian images h are finite; the one for pairing is unsurprisingly the 
OAI prescription, while that for SO (7) is the opposite, orthogonalizing states of low seniority 
against states of higher seniority. These finite, Hermitian images are: 



^pairing _ ^^f^ _ ^t^t^^^ 

{yt + b)sh + S -d- sUhs. 



(80) 
(81) 



Note that these general Hermitian images are not simply related to the truncation of the full 
image through an overall renormalization, as is the Hermitian Dyson image ([Til). For the 
pairing interaction, the coefficient of the s^s^ss term changes sign, while off-diagonal terms 
such as d)d)ss, which exist in the image in the full space, vanish, and in some cases like 
SO (8) other diagonal terms such as d)Sdd can also vanish (depending on the model space). 

If one uses an "inappropriate" transform S it can induce an unneeded and unwanted 
many-body dependence. This principle we illustrate in the SO (8) model, by deriving up to 
second order the Hermitian image of V'^ but using Uqai- 

3/n 



jI^^ = (fi + ^)sU + d^.d-(l + 



1-^ 



d^ ®s 



d®d 



2^5 
T 



d®d 



+ 



d^®d^\^ss] (82) 



plus higher order terms which we drop; this is equivalent to the standard OAI procedure 



computed in the 2-boson space. In figure 1 we display the spectra of (|8lD on the left, which 



is the exact S0(7) result, and (p^ ) on the right, taking = 10 and N = 7. The distortions 
in the right-hand spectrum from the exact result, such as the overall energy shift and the 
large perturbation in the third band, indicates the importance of the missing many-body 
terms. In OAI these many-body terms would appear implicitly in the A^-dependence of 
the coefficients for the two-body terms. From the existence of ( |HT| ) we see that this strong 
OAI A^-dependence is, for this case at least, a needless comphcation. Therefore it is possible 
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that some of the A^-dependence of OAI is induced by their choice of orthogonahzation and 
could be minimized with a different choice. We are currently exploring how to exploit this 
freedom to best effect. 



VI. EFFECTIVE REPRESENTATIONS AND IMAGES 

So far we have considered the mapping of a Hamiltonian in a truncated fermion space to 
a boson space. That truncated fermion space, however, may be inadequate for reproducing 
the spectrum of the full space. Consider a single j-shell fermion space. If the interaction 
is dominated by pairing, then truncation to just J = 0, 2 {s,d) pairs is reasonable 0; but 
for other interactions, particularly the quadrupole-quadrupole interaction, the s, d-space is 
inadequate To rectify these shortcomings one must introduce an effective interaction 



theory for boson mappings. 

An important issue is at what stage to introduce effective operators. For example, 
one could start from an effective fermion Hamiltonian. Because the fermion-pair basis is 
non-orthogonal, however, computation of the effective fermion interaction would be tricky; 
furthermore, the starting Hamiltonian would have a number dependence which would be 
difficult to separate from the number dependence induced by the subsequent boson mapping. 
At the other extreme, Sakamoto and Kishimoto [|18| start from a boson image and use 
perturbation theory to account for excluded states. They do not rigorously derive their 
effective interaction and it is not clear that they properly account for the exchange terms, 
etc., included in the norm operator. 

A better, intermediate approach, is to calculate the boson image in a larger space — say 
sdg — and then renormalize the pure sd interaction so as to account for the effect of the 
(7-boson [Q. Ideally, one should start with a sufficient number of species of bosons so as to 
exactly span the fermion Fock space, and then further truncate the space and renormalize. 
The number of species required will depend on the number of bosons (or pairs) A^, however, 
and we know of no prescription for determining this set of bosons (except for A^ = 1 when 
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all are necessary). 

We present a rigorous and general approach to effective interactions. Following the usual 
Feshbach derivation, we partition the boson Fock space using P to project out the allowed 
space and Q its compliment, with P + Q = 1, = P , = Q and PQ = QP = 0. Then 
the truncated representations are simply 



[Hb]^ = PHbP, 

and 1$)^, = -P 1$). Then the generalized eigenvalue equation in the full space becomes 



eff 



B 



eflf 



T 



(Ex) 1$ 



(83) 
(84) 

(85) 



with 



n 



B 



eff 



(E) = PHbP + PHbQ ^ 

^ ^ ^ ^Q{EArB-HB)Q 



-PHbQ 



E 



Q{ENb - 'Hb)Q 

+pMbQ 



QUbP - PUbQ 



E 



q{eUb - nB)Q 



QHbP 
QHbP 



E' 



Q{ENb-'Hb)Q 



QMbP + EA{E), 



(86) 



■ ^ 1 eff 

Nb]^ {E) = PMbP + A{E). 



One can also in principle construct energy-independent, but non-Hermitian, effective rep- 
resentations ||2H|. There is some ambiguity in the definition of the effective representations 



B 



eff 



to be simply PTIbP with the 



eff 



A/b ^ (-E), and in principle A{E) 



as denoted by A{E). For example, one could define 
remaining terms in (^) absorbed into the definition of 
could be anything at all. 

Now consider boson images, where one divides out the norm operator. We suggest that 
effective operator theory may be more efficient when applied to representations rather than 
images, by which we mean that the corrections are smaller. Suppose one started with the 
image Hamiltonian in the full space, Hb as defined previously, and from that constructed 
an effective image in the usual way. 
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Hf = PHbP + PHbQ-^^^^Jj-^QHbP = [Hslr + ^Hb{E). (87) 

Now compare that with the effective image constructed from effective representations, 

hf{E) ^ ([a^^];^ {E))'" [HbT^ ([aT^];^ {E))'" ^hB + AhB{E) (88) 

(leaving aside the issue of the choice of an overall unitary transformation U). Now fiB 
[HB]rp. Which approach is better? In those cases, such as S0(8) and Sp(6), where, the 
P-space decouples completely from the Q space, AHb = but AHb cannot be zero. Hence 
the corrections AHb from using effective representations can be smaller than the corrections 
AHb determined from performing effective operator theory directly on the image. 

We now would like to speculate on the possible use of the ambiguity operator A{E) . In 
effective operator theory the eigenstates are no longer orthogonal because of the truncation 
of the model space; this is expressed by the fact that one uses either energy-dependent 
or non-Hermitian effective interactions. We propose that this non-orthogonality could also 
be embedded in the choice of A{E), so that the similarity transform on the basis is now 
A/b ^ + A{E) |$a). The ambiguity operator A{E) could be chosen so as to minimize 
the energy dependence of the final boson image. Although the similarity transformation is 
now energy-dependent, this would not show up in the calculation of the spectrum, but only 
in the calculation of effective transition operators. These speculations need to be explored 
in greater detail. 



VII. SUMMARY 

In order to investigate rigorous foundations for the phenomenological Interacting Boson 
Model, we have presented a rigorous microscopic mapping of fermion pairs to bosons, paying 
special attention to exact mapping of matrix elements, Hermiticity, truncation of the model 
space, and many-body terms. First we presented new, general and compact forms for boson 
representations that preserve fermion matrix elements. We then considered the boson image 
Hamiltonian which results from "dividing out" the norm from the representation; in the full 
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boson Fock space the image is always finite and Hermitian; in addition we discussed several 
analytic cases for truncated spaces where the image is also finite and Hermitian. Next, we 
gave a prescription which is a generalization of both the OAI and democratic mappings; 
in the most general case for truncated spaces the Hermitian image Hamiltonian may not 
be finite but we have demonstrated there is some freedom in the mapping that one could 
possibly exploit to minimize the many-body terms. This freedom, which manifests itself 
in a similarity transformation that orders the orthogonalization of the underlying fermion 
basis, depends on the Hamiltonian. Finally, we discussed effective operator theory for boson 
mappings. 

This research was supported by the U.S. Department of Energy. The boson calculations 



for Figure 1 were performed using the PHINT package of Scholten [Id 
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FIGURES 

FIG. 1. Spectrum of S0(7) interaction, for 7 bosons, in S0(8) model with exact (left) and 
approximate (right) two-body boson Hamiltonians. 
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